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ABSTRACT 



Aims. A new method to constrain the cosmological equation of state is proposed by using combined samples of gamma-ray bursts 
(GRBs) and supernovae (SNela). 

Methods. The Chevallier-Polarski-Linder parameterization is adopted for the equation of state in order to find out a realistic approach 
to achieve the deceleration/acceleration transition phase of dark energy models. 

Results. We find that GRBs, calibrated by SNela, could be good distance indicators capable of discriminating between cosmological 
models and ACDM model at high redshift. 
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1. Introduction 

^ From an observational viewpoint, one of the fundamental goals 
' of cosmology is to measure cosmological distances and then to 
build up a suitable and reliable cosmic distance ladder. This is- 
sue has recently become even more important due to the evident 
I degeneracy of several dark energy models with ACDM, desp ite 
the advent of the so-called Precision cosmolosv, dEUis 1999h . 

In the last two decades, a class of accurate standard can- 
dles, the supernovae type la (SNela), has been highly studied 
and the results obtained from the use of these objects led to the 
surprising discovery of the acceleration of the cosmic Hubble 
' flow (for a review see Kowalski et al 2008). However these ob- 
jects are poorly detectable at redshifts higher than ~ 1.5, so we 
need distance indicators at higher redshifts in order to remove 
the degeneration of dark energy models aff'ecting current cos- 
mological models (ACDM is a good approximation of the ob- 
served Universe, even though there is still no theoretical basis 
about the nature of its components, but the issue of global evo- 
lution is far from being addressed; for a comprehensive review 
see Copeland et al 2006). A possible solution could be found by 
adopting gamma ray bursts (GRBs) as distance indicators. 

GRBs are the most powerful explosions in the Universe: the 
most likely scenarios for their generation are the formation of 
massive black holes or the coalescence of binary stellar sys- 
tems. These events are observed at considerable distances, so 
there have been several eff'orts to frame them into the standard of 
the cosmological distance ladder. In the literat ure, there ar e sev- 
eral models that account for GRB formation (iMeszaros 2 006). 
All these scenarios involve a similar shock phenomenon: a "fire- 
ball", possibly be supported by a further jet emission. However 
none of these models is intrinsically capable of integrating all 
the observable quantities. 



Despite the poor knowledge of the GRB mechanism, it 
seems that GRBs could be used as reliable distance indicators. 
There exist several observational correlations among the photo- 
metri c and spectral properties of GRBs to support this possibil - 
ity (Basilakos & Perivolaropoulos 20081; iGhirlanda et a l 20061). 
Nevertheless the origin of these spectroscopic and photomet- 
rical correlations is not known very well and there are sev- 
eral eflforts to interpret the behavior of G RB features in a co- 
herent way, by relati vely simple scenarios (iDainotti et al 20081: 
"Ghisellin i et al 200 8). Succeeding in explain the mechanism that 
generates GRBs is one of the objectives of modern astrophysics 
and to clarify these observed correlations in this context would 
make GRBs reliable distance indicators. A complete review of 
the existing luminosity relations for GRBs can be found in 
Schaefer (2007). 

In this paper, we consi der two relations, the one by Liang- 
Zhang (LZ) (Liang & Zhan g 2005b . and the one by Ghirlanda 
(GGL) ( Ghirlanda et al 2004). They are the only 3 -parameter re- 
lations known and have less scatter with respect to the theoretical 
best fit than the other 2-parameter ones. However calibration of 
the relations used has been necessary in order to avoid the cir- 
cularity problem. This means that all the relations need to be 
calibrated for each set of cosmological parameters. Indeed, all 
GRB distances, obtained in a photometric way, are strictly de- 
pendent on the cosmological parameters since, currently, there is 
no low-redshift (z up to 0.2-0.3) set of GRBs available to achieve 
a cosmology-independent calibration. In order to overcome this 
difficulty, Liang et al. (2008), proposed a method in which sev- 
eral GRB relations have been calibrated by SNela. Supposing 
that these relations work at all redshifts and that, at the same 
redshift, GRBs and SNela have the same luminosity distance, it 
becomes possible, in principle, to calibrate the GRB relations us- 
ing an interpolation algorithm. In this way, it becomes possible 
to build a GRB -Hubble diagram by calculating the luminosity 
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distance for each GRB with the well-known relation between the 
luminosity distance di and the energy-flux ratio of the distance 
indicators. 

In the literature there are several paper that use similar 
methods to constraint the cosmological parameters of the 
Concordance Model using GRBs as extension of the SNela 
Hubble Diagrani , (Firmani et al. 2005; Firmani et al. 20061: 
I Wang & Dai 20061: IWang et al. 2007; Li et al. 2008; Li et al.1) . 

Here we take into account a cosmological EoS working 
at any redshift, using GRBs as tracers and adopting again the 
Chevallier-Polarski-Linder (CPL) parameterization. In particu- 
lar we discuss a method which should allow us to obtain an ana- 
lytic cosmology-independent formulation of the luminosity dis- 
tance and then of the distance modulus. After a brief introduc- 
tion to the GRB luminosity relations, we show fits of the data 
obtained by these relations and the results and perspectives of 
the approach are discussed in the last Section. 



2. The theoretical framework 

Our goal is to obtain an analytic formulation of the Hubble dia- 
gram valid at any redshift. We start from the Friedmann equation 



SttG 



kc^ 



(1) 



We obtain, by some algebra, the following equation in terms of 
the density parameter 



Qo(^P-(no-l)(^y 



(2) 



where the subscript indicates the present value of the pa- 
rameters. From now onwards, we take into account a spatially 
quasi-flat Universe, k ^ 0; the contribution of the curvature 
will be negligible and we have Do ~ 1, as suggested by the 
latest CMBR ( Komatsu et al 2008) and the SNela observations 
(iKowalski et al 200'8|) . However, in the final section, we will per- 
form a test to verify this assumption with observations coming 
from GRBs. Now if we translate in terms of redshift z. 



— = 1 +z, 
a 



the previous equation reduces to 



3(w+l) 



(3) 



(4) 



The w-parameter indicates the EoS w = p/p, where p and p 
are the pressure and the matter-energy density of the Universe, 
respectively. Considering the CPL parameterization of the EoS, 
(pievallieretaL2001): 



w(z) = Wo + w« 



and substituting into Eq.©, we obtain: 



Hiz) = Ho 



Ud+z) 



exp 



(5) 



(6) 



which enters directly in the expression of the distance modulus 
^i{z) = -5 + 5logd,{z), (7) 



where di{z) = c(l + z)Di{z) and where 
Diiz) 



Jo mo' 



(8) 



This means that an analytic expression for n can be achieved. 
The integral Di in Eq.© can be solved giving a Gamma function 
of the first kindQ: 



exp | — IT 



1 -h 3wo + ?>Wa 



1 +3^0 +3^0 
2 

^=0 



(9) 



Substituting such an expression in the distance modulus, we 
obtain a model for data fitting which could work, in principle, 
at any z- The obtained expression for the Hubble parameter H{z) 
is independent of the density parameters, Dm and Da, provided 
that their sum is equal to 1 . 

We use the CPL parameterization not only for the dark en- 
ergy component, but for the total energy-matter density of the 
Universe. This assumption works because dark and baryonic 
matter contribute with a null pressure while the radiation com- 
ponent is negligible in matter- and dark energy-dominated eras. 
Furthermore, the analytical formulation that we adopt for the lu- 
minosity distance is assumed valid at any redshift z. 

3. GRB luminosity relations 

In the last years, thanks to several spacecraft missions capa- 
ble of observing this high energy region, the main features of 
GRBs have become better known. Recently, some photomet- 
ric and spectroscopic relations between GRB observables have 
been found and then the hypothesis that these objects could 
be considered suitable distance indicators has become feasible. 
Nevertheless, there is no theoretical model that fully explains 
these relations so the GRBs cannot be considered as standard 
candles. For a detailed review of the observational features see 
Schaefer (2007). 

Here, we take into account the existing 3 -parameter re- 
lations. This choice has been made because these rela- 
tions place better constraints on the data giving less scat- 
ter between th e theoretical relation and the experimental data 
(IScha efer 2007). The first relation is the so-called Liang-Zhang 
relation, ( Liang & Zhang 2005), which allows us to connect the 
GRB peak energy, Ep, with the isotropic energy released in the 
burst, Eiso, and with the jet break - time of the afterglow optical 
light curve in the rest frame, measured in days, t^, that is 



logEiso = a + bi\og 



Ep{\+z) 



+ b2 log ■ 



tb 



(10) 



?>mkeV ^(l-hz)lJaj 

where a and hi, with / = 1,2, are calibration constants. 

The other re lation is that given by Ghirlanda et al 
( Ghirlandaet al 2004). It connects the peak energy Ep with the 
collimation-corrected energy, or the energy release of a GRB jet, 
Ey, where Ey = FteamEiso = 1 -cos(e)Eiso, with 6 jet the jet open- 
ing angle defined in Sari et al (1999): 



Ojet =0.163 



mrjy 



3/8 

Z/ \Eiso,52 



1/8 



(11) 



^ In our case, the variable of the Gamma function, z, is always posi- 
tive so that we have no problem of discontinuity in applying the gamma 
function in the following calculations. 
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Table 1. Results of the fits. SNela is only for the supernovae 
data, LZ is for the GRBs data obtained from the Liang-Zhang 
relation, GGL for the Ghirlanda et al. one. 



Relation 


Wq 










SNela 


-0.910 


+ 


0.070 


0.755 


+ 


0.054 


0.983 


LZ 


-1.39 


+ 


0.38 


1.18 


+ 


0.37 


0.817 


GGL 


-1.46 


+ 


0.38 


1.36 


+ 


0.32 


0.812 


LZ + SNela 


-1.15 


+ 


0.10 


0.93 


+ 


0.11 


0.933 


GGL + SNela 


-1.42 


+ 


0.12 


1.24 


+ 


0.13 


0.920 



where Eiso,52 = Eisol^^ ergs, is the circumburst particle 
density in 1 cm"^, and riy the radiative efficiency. The Ghirlanda 
et al. relation is 



log Ey = a-\- b log ^ — , 



(12) 



where a and b are two calibration constants. 

From these relations, we can directly obtain the luminosity 
distance di from the well-known formula which connects di with 
the isotropic energy Eisq and the bolometric fluence S ijoio '• 



(13) 



bolo 



from which it is easy to compute, for each GRB, the distance 
modulus [I = and its error given by (Liang et al 2008.) : 



{2.5criogEi,af + '^^^CTS tola bolo f 



(14) 



with o-\ogEi,„ and cr^^^/^ obtained from the error propagation ap- 
plied to Eq.JTOl) and Eq.([T2l). Moreover, we assume that the error 
in the determination of the redshift z is negligible, as well as for 
the radiative efficiency rjy. We note also that the assumption of 
a well-known hq is a strong hypothesis since the goodness of 
the fits depends, in particular, on this parameter. The GRB data 
sample is taken from the already cited work by Schaefer. We take 
into account 27 events with extremely precise data. This sample 
is the same one adopted in Capozziello & Izzo (2008). 




Fig. 1. Redshift-distance modulus diagram for the GRB-^SNela 
sample. The black dots are the GRBs , the blue ones are the 
SNela. The red line is the best fit obtained from the data, with the 
dashed line representing the confidence limits at 3cr. The error 
bars on the supernova data are not represented because they are 
negligible. 




Fig. 2. Comparison between the best fit of /j, and the observed 
distance modulus fiobs at any redshift. The black dots are the 
GRBs data and the red line is the best fit curve representing the 
theoretical distance modulus. 



4, The data fitting 

The next step is the fit of the GRB sample with the empiri- 
cal relations, Eqs.([TQb,([T2b, described in Sect. 3. The aim is to 
achieve an estimate of the CPL parameters and consequently to 
determine the trend of the EoS at any redshift, using the an- 
alytical relation, Eq.©. We are considering the same sample 
of 27 GRBs used in Capozziello & Izzo (2008) in which we 
have added the sampl e of SNela by the Union Supernova Survey 
(iKowalski et al 2008). 

The numerical results of the fits are shown in Table [TJ where 
we obtain a robust estimation of the CPL parameters for both 
the relations used, with and without SNela data. An immediate 
comparison is done with the best fit applied only to the SNela 
sample. It is evident how adding GRBs to SNela data completes 
the knowledge of and the accuracy on the EoS parameter w. 



In order to measure the goodness of the fit, we use the test, 
see Table [T] The R^ test is a measure of how successful the fit is 
in explaining the variation of the data (see for details Draper & 
Smith 1998). An R^ close to 1 .0 indicates that we have accounted 



for almost all of the variability with the data specified in the 
model. As a standard, the R^ test is the square of the correlation 
between the response values and the predicted response values, 
that is: 



R' 



1 



SSE 

Jsf 



1 



(15) 



where SSE is the sum of the squares due to errors and it mea- 
sures the total deviation of the response values from the fit and 
S ST is the sum of squares about the mean: y is the predicted 
response value, y is the mean value and the w/ are the weights on 
the values. 

The extension of the supernova Hubble Diagram with the 
GRB data can be used to improve our knowledge of the trend at 
high redshift. In this way, also using the GRB data, we show in 
Fig. [3l the distance modulus ji versus the redshift z, in a log- 
arithmic scale. The best fit curve, obtained with Eq.®, is also 
reported. A more detailed analysis confirms the presence of a 
transition (re- acceleration) redshift around z = 0.5. 
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Redshiftz 

Fig. 3. Redshift-distance modulus diagram for the GRB+SNela 
sample versus redshift in logarithmic scale. 



In Fig. [21 we plot the comparison between the theoretical 
[ith and the observed distance modulus [lobs at any redshift, the 
residual plot. A smooth trend up to z ~ 3.5 in the residual curve 
can be immediately detected. Beyond this limit, we have 3 GRBs 
that exceed, by the same side, the 3(T confidence limit of the best 
fit. This discrepancy is clear in Fig. [3l where we plot the best fit 
for the combined sample in the case of an LZ relation with a 
logarithmic scale for the redshift. 

This fact is fundamental for the goodness of the fit because 
these GRBs represent the most distant objects that one can use to 
make such an analysis and their weight on the fit is very high, in 
the sense that they appear to be not accurate distance indicators. 

However there is strong evidence that the data of these 3 
GRBs, GRB 050505, GRB 050904 and GRB 060210, reported 
in the Schaefer catalog, are uncertain. For the first and the la st 
GRB the peak energies are underestimated (ICabrera et al. 20071) . 
As a consequence we would obtain an underestimated value for 
the bolometric fluence. GRB 050904 is the most distant GRB 
considered and it shows some probelms related to the peak en- 
ergy reported by diff'erent authors, that diff'er by a factor of 
~ 3. Moreover the afterglow of GRB is complicated by flares 
and re-brightenings so that the standard afterglow model gives 
an extremely h igh value of the circumburst medium density 
n ~ 700cm"^ (iFrail et al. 2006h . contrary to the assumed value 
n = 3cm"^. 

For this reason we repeat the analysis described above with- 
out these 3 GRBs, obtaining a better value than the previous one 
for the test. The results of these corrected fits are shown in 
Table [3] In Fig. IH we plot the best fit with this corrected sam- 
ple. From these results we conclude that the complete sample 
gives diff'erent results than does corrected sample, the first one 
suggesting a phantom/quintessence regime for the present epoch 
while the second one fits an accelerating ACDM model. This 
last result is confirmed by the following analysis, where we have 
performed a Monte-Carlo-like procedure for the comparison of 
the results with the usual likelihood estimator given by 



X 



(j^th(Zi) - flobs(Zi)f 



(16) 




1.5 2 
Redshift z 



Fig. 4. Redshift-distance modulus diagram for the corrected 
GRB-hSNela sample. The red line is the best fit obtained from 
the data, with the dashed line representing the confidence limit 
at 3o-. 

modulus uncertainty. The results of this analysis are shown in 
Table[2l where we can see the improvement obtained by the GRB 
sample corrected for the 3 "wrong" GRBs. 

We have adopted a similar procedure in the case of an EoS 
evolving with redshift and where jith is obtained by Eq.®. The 
result of this analysis is plotted in Fig. [5]where the best fit value, 
the cross in the figure, corresponds to the value wq = -0.84+0. 14 
and Wa = 0.72 + 0.06, and where the boundaries correspond 
to Icr, 2cr and 3cr confidence levels, in a good agreement with 
the results obtained, see Table [3j using our theoretical relation, 
Eq.©. 

From this analysis, we conclude that the corrected sample 
agrees fairly well with the ACDM model with a small contri- 
bution of the curvature parameter, k = 0.01 + 0.04. Thus, the 
method delineated in Sect.2 seems a good approximation of the 
observed cosmography and agrees very well with the ACDM 
model, so that we can argue that GRBs could be good distance 
indicators at redshift values up to z = 4. 



in the context of a ACDM model of the Universe, where /dth is 
the distance modulus computed from the Eq.(|7]) and Eq.©, Zi is 
the observed redshift for each GRB and ct/ the observed distance 



5. Discussion and conclusions 

Starting from the Friedmann equation, we have investigated a 
new method to constrain the cosmological equation of state at 
high redshifts. In particular we obtain an analytical formula for 
the distance modulus so we could directly estimate the parame- 
ters of the cosmological model considered. The working hypoth- 
esis involves the use of GRBs as distance indicators at high red- 
shift, well beyond the distance where SNela have been detected 
to date. The CPL parameterization for the EoS has been explic- 
itly used for the whole matter-energy content of the Universe as 
a suitable approach to investigate the parameter w = w(z) and 
discriminate values with respect to the ACDM model. In partic- 
ular, regarding the Friedmann equations, we have obtained, in 
the case of the LZ relation, the epoch for the transition between 
the deceleration-acceleration phases at a redshift value of z ~ 5, 
with a reliable confidence level. This is a value that, if higher 
than the redshift of the farthest GRB used, could be in agree- 
ment with current quasar formation scenarios. Also, we are in 
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Table 2. Cosmological density parameters, with uncertainties computed at Icr confidence Hmit, obtained by a Monte Carlo procedure 



Sample 










x' 


UNION H 


hGRB 


0.26 + 0.14 


0.73 + 0.14 


0.01 + 0.04 


1.032 


UNION H 


h ORB corrected 


0.25 + 0.10 


0.74 + 0.135 


0.01 + 0.035 


1.00027 




-1.4 -1.3 -1.2 -1.1 -l.O -O.S -0.7 -0.6 



Fig. 5. 68%, 95% and 98% constraints on wo and w« obtained 
from the UNION sample and the GRB sample corrected for the 
3 outlying GRBs. The cross represents the best fit value and it is 
in a good agreement with that found using the theoretical model 
described in Sect.2. 



good agreement with the observed phantom/quintessence regime 
at the present epoch, that is for z ^ 0, we obtain w < - 1 . 

So we reject the current phantom regime by this analysis, ob- 
taining for Wo a value in agreement with the ACDM model at the 
present epoch. The method, while preliminary, seems to indicate 
that GRBs could be used as standard candles once a reliable uni- 
fied model of their photometric and spectroscopic quantities is 
achieved (some relevant results are presented in Ghisellini et al 
2008). However, more robust samples of data are needed and a 
more realistic EoS (with respect to the simple perfect fluid mod- 
els) should be taken into account in order to suitably track red- 
shift at any epoch (see for example Capozziello et al 2006). 

With improving observations, in particular with the launch of 
new satellites devoted to GRB surveys, such as Fermi- glastEI 
and AGIL^ one should be able to expand the samples of GRBs, 
possibly with data coming from objects at higher redshift. 

Considering these preliminary results, it seems that GRBs 
could be considered as a useful tool to remove degeneracy and 
constrain self-consistent cosmological models. The matching 
with other distance indicators would improve the consistency of 
the Hubble distance-redshift diagram by extending it up to red- 
shift 6-7 and higher. 



Table 3. Results of the fits corrected for the 3 "outying" GRBs. 
SNela is for the supernova la data, LZ is for the GRBs data ob- 
tained from the Liang-Zhang relation, GGL for the Ghirlanda et 
al. one. 



Relation 








LZ + SNela 


-0.95 + 0.01 


0.74 + 0.01 


0.999 


GGL + SNela 


-0.865 + 0.005 


0.66 + 0.005 


0.999 
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